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Abstract 

We present a geometric, model- independent, argument that aims to explain why the Tsal- 
lis entropy describes systems exhibiting "weak chaos" , namely systems whose underlying 
dynamics has vanishing largest Lyapunov exponent. Our argument relies on properties of 
a deformation map of the reals induced by the Tsallis entropy, and its conclusion agrees 
with all currently known results. 
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1. Introduction 



The Tsallis entropy is a functional first introduced in the Physics hterature in 1988 
[1], providing an alternative to the Boltzmann/Gibbs/Shannon (BGS) entropy. It aims to 
describe collective phenomena with long-range spatial and temporal correlations, systems 
whose phase portraits exhibit a fractal-like behavior for which there is no reason why the 
description through the BGS entropy should be accurate or even valid [2] (and references 
therein). Although the dynamical basis of the Tsallis entropy remains unclear so far, 
not totally unlike that of the BGS entropy [2]- [6], there is certain progress that has been 
made in identifying characteristics of systems that seem to be well described by it. One 
class of such systems are ones exhibiting "weak chaos" or being at the "edge of chaos" for 
some subset of their parameter space, more concretely referring to systems with vanishing 
largest Lyapunov exponent [2], [7]- [10]. 

A general and relatively formal justification of why such systems are described by the 
Tsallis entropy has been lacking so far [2], [11]. This is the issue that we are addressing in 
the present work. To be more specific, we present an argument about the converse: if a 
system is described by the Tsallis entropy then its largest Lyapunov exponent vanishes. 
This fact is in agreement with all currently known results pertaining to the microscopic 
dynamics of systems described by the Tsallis entropy. As a by-product we also comment 
on modified definitions of asymptotic numerical invariants of dynamical systems better 
reflecting the properties of the Tsallis entropy [2], [11]-[16]. Our approach relies on the 
metric aspects of Riemannian manifolds.. For this reason some arguments can also be 
applied, or non-trivially generalized, in the much more general classes of (geodesic) metric 
spaces, such as spaces with a (negative) upper bound A; < on their curvature, known 
as CAT{k) spaces. A feature of our presentation that sets it apart from most other 
descriptions of Riemannian methods addressing a Physics audience, is that it stresses, 
as much as possible the synthetic, rather than the analytical viewpoint. The advantage 
of this approach is that it is immediately generalizable to more general geometric struc- 
tures such as the CAT{k) spaces, mentioned in this work, that do not have very "nice" 
smoothness/regularity properties. In Section 2, we provide a very brief background in Rie- 
mannian geometry in an attempt to make the exposition at least minimally self-contained 
for our audience. We forego essentially all details and infinitesimal calculations as they 
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can be found in any exposition of Riemannian geometry, such as [17], [18]. In Section 3 
we present our argument. Section 4 points toward possible extensions of the geometric 
structures and future research. 

2. Geometric preliminaries 

Riemannian manifolds M of dimension n are structures that are made up of lo- 
cal patches of Euclidean spaces of (integer Hausdorff) dimension n. Each of these 
patches can be parametrized by a coordinate system. The patches are "glued" together 
via functions having desired regularity properties (homeomorphisms, PL-equivalences, dif- 
feomorphisms etc). For most physical applications Riemannian manifolds are assumed to 
be smooth enough, namely the transition functions that glue together local patches of R" 
are assumed to be C°°. Riemannian manifolds are therefore metric spaces whose metric 
is locally Euclidean since it is descended from the Euclidean metric of each patch. Quite 
importantly, they are length spaces and even more specifically, they are geodesic spaces 
(also called intrinsic length spaces), namely spaces for which the distance between any 
two of their points is realized as the infimum of the lengths of all piecewise smooth curves 
joining them. A direct consequence is that Riemannian manifolds have built-in the con- 
cept of direction, as can be seen by the fact that their metric structure is determined by 
specifying an ellipsoid in the neighborhood of each of their points. This ellipsoiid is a re- 
minder of the Euclidean patches that were put together to form the metric of Riemannian 
manifolds. An ellipsoid is analytically described as the unit sphere of a positive-definite, 
non- degenerate quadratic form g, called the metric tensor, which is therefore defined 
though its sections with each 2-dimensional subspace of the tangent space T^M at each 
point X G M. The tangent space T^M is the set of vectors that are tangent at a; G M 
to all curves within M passing through x. Hence, roughly speaking, elements of T^M 
are derivations in the space of smooth real functions defined on the neighborhood of each 
X G M. Equivalently, the curves passing through x are the integral curves of the 
corresponding elements of T^M. One can map each element X G T^M to its corre- 
sponding integral curve c{t) by using the exponential map exp^ : T^M M defined 
by exp,(X)(t) = c(t). 

The fact that the two metric structures, one which is induced by the underlying topol- 
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ogy of M" and the other which is induced by the length structure are equivalent, is 
non-trivial, and is the content of the Hopf-Rinow theorem [17], [18]. In that theorem it is 
assumed that the Riemannian manifold is complete, meaning that no points have been 
eliminated from it, as for instance in the case of M^\{0} where denotes the origin. 
This is more compactly expressed by requiring exp^ to be complete. The process of 
actually requiring the exponential map to be complete may appear quite technical and 
far removed from the needs of Physics, but such a statement would only be partly cor- 
rect. There is at a least one notable case of use of non-complete manifolds for physical 
purposes and this the case of the singularity determination in black holes. In the present 
work we will ignore cases of non-complete manifolds since they may be encountered rarely 
in Statistical Mechanics so they are not very important for our purposes. 

It may also be noticed that there is nothing really "sacred" about using a non- 
degenerate quadratic form, namely an ellipsoid, to determine the distance function as 
is done in Euclidean/ Riemannian spaces. One could use an appropriately scaled version 
of any symmetric convex body instead, which is the unit ellipsoid of a general Banach 
norm. This is indeed done in the case of Finsler manifolds. The general determination 
of a symmetric convex body, however, requires infinite parameters, even in the planar 
case. By contrast, the determination of an ellipsoid requires knowledge of just the lengths 
of its axes. Moreover, no other Banach spaces except the square summable/integrable 
cases are induced from an inner product. The inner product is a very well- understood 
and convenient technical device, so we would like to keep using it in our models, unless 
particular requirements force us to abandon it. 

As in any metric spaces, also in the case of Riemannian manifolds M, one would like 
to find ways of determining the distances between any two of their points. This turns out 
to be a hopeless task, even in principle. First of all, to appropriately determine a distance 
d{x,y), x,y & M, the corresponding Riemannian metric g will have to be determined. 
This amounts to determining n{n + l)/2 numbers at every point of an n- dimensional 
Riemannian manifold as g is a symmetric, non-degenerate quadratic form. This is a 
local, even though quite often cumbersome, calculation so it is feasible. Determining, 
however, the minimum distance over all (rectifiable) curves between any two points x 
and y is an intractable task. We can try, for instance, to find the distance between two 
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matrices belonging to the (linear) representation of a Lie group, whose set is a Riemannian 
manifold when endowed with the Cartan-Killing metric. This turns out to be not a very 
simple task. To see how things can get even harder, consider the following example: let 
iV be a Riemanian manifold which is not simply connected, meaning that there are loops 
in N that are not continuously deformable ( "homotopic" ) to points. It turns out that 
each such N has a "universal cover" which is simply-connected. Although globally 
and A^ may be appear different, they are locally (arc- wise) isometric, in the neighbor- 
hood of (almost) every point. Isometry signifies that the distances are locally preserved 
in going from A^ to A^. To be more specific, the length structure lifts from A^ to A^ by 
assigning to each curve of N its image in A'". What is still unclear is how to determine 
the metric on N. Someone may ask, for instance, whether the diameter of N defined 
as sup d{x,y), \/ x,y G N is finite or not. This reduces to deciding whether the funda- 
mental group 7ri(A^), which is an algebraic structure (a group) measuring the algebraic 
complexity of families of homotopic loops in A^, is finite or not. This is however an unde- 
cidable proposition, therefore the diameter of N cannot be effectively calculated [19], [20]. 

To get around these difficulties, one would fike to define effectively computable invari- 
ants of M which are local in the sense that they are defined for a neighborhood of each 
point of M. Such invariants are usually related to the Riemann curvature tensor, which 
is a relatively complicated non-linear expression of the metric tensor g its inverse g"^ 
and their second derivatives [18], [19]. As such an expression is not very illuminating and 
it manifestly depends on the system of coordinates employed, one may wish to specify it 
in an alternative way, which would exhibit its linear-algebraic character and which would 
make more prominent its coordinate independence, and consequently its geometric mean- 
ing. This is a partial motivation for introducing the concept of a Levi-Civita connection 
[17], [18]. A connection, as the word indicates, is a way of splitting the tangent space of 
a Riemannian manifold into a horizontal and a vertical part T^M — (T^MY + {T^M)-^ 
at each point x e M. This is an orthogonal splitting which is required to be compatible 
with the properties of g. The horizontal part of this splitting is used to compare or 
"connect" the tangent spaces over different points of M, hence the name of the struc- 
ture. The connection amounts to essentially choosing how to identify the tangent spaces 
T^M at different points x e M by declaring which curves are horizontal, by locally 
specifying their tangent vectors. Hence, a connection can be also seen as a derivation on 
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TxM, or cquivalcntly as the directional (covariant) derivative of vector fields Y along 
curves with tangents determined by the vector field X. The notation for this is Vx^- 
This (covariant) derivative operation can be extended to tensor products of T^M and 
its duals T*M by linearity and Leibniz's rule. The Riemann curvature tensor expresses 
the path dependence (or lack of integrability) of the transport of a vector around an in- 
finitesimal loop. To determine it, one transports a vector Z G T^M around the loop 
with the shape of an infinitesimal rectangle whose sides are the vector fields X,Y E TM. 
Following this motivation and wanting to keep the geometric character (reparametrization 
independence) of the resulting object by requiring it to be a tensor, the Riemann tensor 
is defined as 

R{X,Y)Z = Vx^yZ -VyVxZ -Vix,Y]Z (1) 

Here, the square brackets indicate a commutator. From this relation one immediately 
sees that the Riemann tensor is the obstruction to the integrability of the vector Z 
across the distribution determined by X and Y and, as such, it quantifies curvature. 
Hence, the curvature is fundamentally a 2-dimensional concept, although it appears far 
more complicated when seen in the totality of all 2-dimensional subspaces of T^M, 
formally known as the Grassmann manifold Gn,2iM). To quantify more concretely 
and geometrically the concept of curvature, several scalar quantities derived from the 
Riemann tensor are utilized [17]- [20]. The most general of the widely used ones, is the 
sectional curvature K, which is defined in each 2-dimensional subspace of T^M. Let 
such a subspace S be spanned by the orthonormal vectors 61,62 with respect to the 
Riemannian metric g. Then 

Kj:^g{R{ei, 62)62,61) (2) 

The geometric interpretation of is provided by the following well-known fact [18], [20]. 
Consider a circle of radius r in E. Use the exponential map exp to map this circle 
to a closed curve of length Ir in M. Then 

3 27rr - Ir , . 

= - lim — (3 

TT i — ^0 r 

We immediately see that positive sectional curvature corresponds to the actual shrinking 
of the length of c^, whereas negative sectional curvature is equivalent to an expansion 
of Ir- This is the critical geometric observation that will be used in the sequel and on 
which our whole argument relies. Conversely, if one knows the sectional curvature along 
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all 2-dimensional planes at a point x G M, then one can determine the Riemann tensor by 
the algebraic operation of polarization [17], [18]. It may be worth mentioning at this point 
that the definitions of the Riemann tensor through the connection or through a metric are 
equivalent. This happens because Riemannian manifolds have the following nice property: 
there is a unique symmetric connection, which is compatible with the metric g, meaning 
that it preserves the inner products of any two vectors X,Y E T^M. As a reminder, 
a connection is symmetric, also called torsion-free, if it satisfies VxY — VyX = [X, F], 
so in a coordinate basis the corresponding connection coefficients (Christoffel symbols) 
are symmetric. For completeness, we state the explicit expression of the connection in 
terms of the metric g, whose construction was inspired by group cohomology, and is due 
to Koszul 

2g{VxY,Z) = X[g{Y,Z)]-Z[g{X,Y)]+Y[g{X,Z)]-g{X, [Y,Z])+g{Z, [X,Y])+g{Y, [Z,X]) 

(4) 

In Riemannian geometry a curve is called geodesic between two points x,y E M, if it 
is a stationary point of the length functional on all (rectifiable) curves joining x and 
y. The concept of geodesic reflects the Euchdean concept of the "most straight line" 
between the two points. This line does not have to be a shortest one though, or even to 
be unique. Consider, for instance, a major circle joining two points of a 2-sphere. It is a 
geodesic even though it may not be minimal. Moreover, if these two points are the south 
and the north poles then this geodesic is only one member of an infinite continuous family. 
Prom its variational characterization, or as a definition, it turns out that a geodesic whose 
tangent is X e TM can be analytically determined as the solution of 

VxX = (5) 

which states that the tangent vector X to the curve remains parallel to itself, hence the 
corresponding integral curve is the "straightest possible" line joining x and y [17], [18]. 
In more physical terms, the equation defining geodesies is the Euler-Lagrange equation 
of the kinetic energy of a free point particle moving on M. Therefore the geodesies of 
M are the trajectories of freely moving particles on M. We have tacitly assumed in the 
preceding comments that all geodesies are parametrized by arc-length. 

Another fundamental concept in Riemannian geometry is that of Jacobi fields. The 
idea is that these are vector fields joining nearby geodesies. As such, Jacobi fields describe 
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how nearby geodesies converge or diverge with respect to each other. They are the 
means by which gravitational interactions in General Relativity are determined (albeit 
on manifolds with a "metric" of indefinite/Lorentzian signature). They quantify the 
tidal forces exerted on an object by gravity. So, the Jacobi fields can be variationally 
characterized as extremizing the equation of geodesies. In more physical terms, this 
amounts to setting the second derivative of the kinetic energy of a free particle to zero. 
To find such an extremum one differentiates the geodesic equation using variations of via 
a family of interpolating geodesies. Then the equation for Jacobi fields is derived, known 
as the Jacobi or geodesic deviation equation. It is 

Vx^xJ + R{J,X)X ^0 (6) 

where X is tangent to the geodesic, and where J indicates a Jacobi field [17]. [18]. We 
see that the Jacobi equation is a linearization of the equation for geodesies. It is extremely 
difficult to explicitly solve, in practice, either the geodesic or the Jacobi equation, except in 
few cases of particularly simple Riemannian manifolds. Such a class of simple manifolds, 
which are used as prototypes/models for comparison with most other cases, are provided 
by the manifolds of constant sectional curvature. To keep the topological complications 
to a minimum, it is the simply connected ones (universal covers of all the others) that are 
used as such models. They are called space forms. The idea of comparison geometry is 
the same as in physics: one tries to understand some simple model cases and then uses 
the conclusions reached from them to compare the results with more general spaces or 
realistic systems. Space forms are simple enough that one can actually explicitly calculate 
their sectional curvatures K [17], [18] and find 

R{X, Y)Z = K{g{Y, Z)X - g{Z, X)Y) (7) 

Then following the Cartan-Ambrose-Hicks theorem one can prove that two such space 
forms having equal sectional curvatures are isometric [17], [18]. This is a very concrete 
example that shows that curvature indeed determines the metric properties of a manifold 
which, after all, was a motivation behind the introduction of the Riemann curvature tensor 
[19], [20]. We see that for a manifold with constant negative sectional curvature K < 0, 
the Jacobi equation has the general solution 

n 

J{t) = ^|a,cosh(v^^ t) + biSmh{^/^ t)\ ei{t) (8) 
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where a^, bi are constants, {ei{t)} arc parallel orthonormal vectors (Fermi basis) and 
t is the parameter of the geodesic whose tangent is X{t). We observe therefore that for 
a manifold with constant negative curvatTire, nearby geodesies separate from each other 
exponentially, asymptotically. By comparison, if the manifold is flat namely if the Rie- 
mann tensor/sectional curvature vanishes identically, then the separation between nearby 
geodesies is linear as can be seen again from the Jacobi equation. It is also remarkable 
that there is nothing in between these behaviors, for large values of t. Geodesies of 
Riemannian manifolds either separate linearly {K = 0) or exponentially {K < 0). No 
other asymptotic rate of separation of nearby geodesies is possible, such as a polynomial 
(power-law) separation with exponent different from one, a fact that is obvious from (8). 

3. Zero Lyapunov exponents and Tsallis entropy 

Ideally, one would like to be able to derive the Boltzmann/Gibbs/Shannon form of 
entropy from the dynamics of the underlying system. This is fundamental in Boltzmann's 
approach to systems with many degrees of freedom [3]- [6], although it should be noted 
that such viewpoint is not necessarily shared by Gibbs [21]-[23]. It is clear that this goal of 
Boltzmann has not been realized yet [2], [5], [6]. The process, however, has given a consid- 
erable impetus to the development of measure-theoretical aspects of dynamical systems 
such as ergodic theory [24], [25]. An important realization is that the description of systems 
whose thermodynamic behavior is described the BGS entropy depends on whether they 
possess any (evolution) invariant measures, known as ensembles, in phase space. Such 
measures, whose building blocks are called ergodic measures, provide the stationary prob- 
ability distributions of equilibrium statistical mechanics. The basic idea, due to Gibbs, is 
that the classical ensembles give equivalent thermodynamic descriptions [21]. This how- 
ever turns out not to be true [26]-[28] generally, and especially if the systems possess 
long-range interactions [26], [27]. A broad set of examples of systems possessing ergodic 
measures are provided by the uniformly hyperbolic (Anosov) systems [24], [25]. Uniformly 
hyperbolic dynamical systems are ones for which there is a decomposition of the tangent 
space at each point, varying continuously but not necessarily smoothly, to two transver- 
sal distributions of exponentially distance-expanding and distance-contracting subspaces 
of the tangent space. To be more concrete, a flow [24], [25], is a one-parameter family 
of diffeomorphisms of the Riemannian manifold ("phase space") {M,g) indicated by 
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ft'.M^ M. Let the vector field generating tliis flow be indicated by 

m - j^{Mxm=o (9) 

This flow is called Anosov, if V x G M there exists a direct-sum decomposition T^M = 
E^{x) ® E^{x) © E''{x) and constants c > and < A < 1 < such that Vx e M 
the following four properties hold 

1. E^{x) is a one-dimensional subspace generated by X{t) 

2. E^{x) and E'^{x) remain invariant under the flow, namely dx{ftE^{x)) — E^{ft{x)) 
and d,{ftE-{x)) = 

3. All Y e E^{x) are uniformly contracting, namely < cA*||y||, for t > 

4. All Z e are uniformly expanding, namely < c//~*||Z||, for 
t > 

To determine the exponential rate of expansion and contraction of distances on E'^{x) and 
E^{x) aptly named unstable and stable spaces, of an Anosov system, several quantities 
have been introduced. Probably the simplest of them is the set of the system's Lyapunov 
exponents. This set expresses the asymptotic instability of the evolution of the system 
under inflnitesimal perturbations. In other words, the Lyapunov exponents quantify how 
fast the inflnitesimal perturbations grow or contract in time. Let Y G TM and let its 
norm with respect to the Riemannian metric g be = ^ g(Y, Y). The Lyapunov 

exponent of the perturbation in the direction of Y at x G M along the evolution 
(trajectory) of the system determined by ft is deflned by 

X.{Y)^}im I log \\d,{ftY)\\ (10) 

Prom their definition, it turns out that Anosov fiows having at least one positive Lyapunov 
exponent, exhibit exponential sensitivity to the initial conditions, which is one definition 
of (strong) chaotic behavior [2], [24], [25]. When this behavior is combined with the recur- 
rence properties on their (assumed, for simplicity to be compact) phase space, it results 
in ergodic behavior which ultimately provides the known, though not yet completely 
convincing, dynamical reasons for use of the BGS entropy in describing the system's sta- 
tionary, or nearly stationary, behavior [24], [25]. The way the BGS entropy appears in 
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this context is as the t — )■ oo hmit of another characteristic of the system known as its 
Kolmogorov- Sinai (KS) entropy [24], [25]. The KS entropy is defined in a way similar to 
that of the BGS entropy, after which it is modelled, where in addition the information 
theoretic ancestry of it also becomes apparent [24], [25]. The field of hyperbolic dynamics 
is replete with important technicalities, as such fiows generically result in non-smoothly 
varying continuous unstable and stable distributions, which are the sub-bundles of TM 
generated from E^{x) and E''^{x) through ft. We will however ignore such issues in 
the sequel and we will pretend instead, that the definitions that we use have the required 
degree of smoothness and refer to standard references such as [24], [25] for most of the 
non-trivial regularity aspects. Various generalizations of Anosov flows exist, for instance 
when the norm bounds are not uniform or when the center space, generalizing E^[x), is 
multidimensional and when within which sub-exponential expansions and contractions are 
allowed [24], [25]. However the fundamental ideas behind the definition of Anosov flows 
stated above remain essentially unchanged. 

The most common, arguably, set of examples of such uniformly hyperbolic (Anosov) 
systems is the geodesic flow in the tangent bundle of a manifold of negative sectional cur- 
vature [29] . To see that this Anosov flow indeed possesses at least one positive Lyapunov 
exponent is easiest in the case that (M, g) has constant negative curvature. Then the 
desired result follows straightforwardly from (8) when i — )■ oo: since nearby geodesies 
deviate from each other exponentially in the direction of a vector Y, the corresponding 
Lyapunov exponent which quantifles this exponential deviation, according to (8) must 
be positive. Proving this for the geodesic flow on a manifold of non-constant negative 
sectional curvature is non-trivial though. Such a statement can be made plausible by 
employing the Rauch comparison theorem [17]: it, roughly, states that for two Rieman- 
nian manifolds Mi and M2 of not necessarily the same dimension, if the sectional 
curvatures along all corresponding 2-dimensional planes Ei and E2 at each correspond- 
ing set of points obey K^-^ > Ky,2, then the corresponding Jacobi fields Ji and J2 
which are normal to the respective radial geodesies, which are assumed to be arc-length 
parametrized, obey [1^2(^)11 > ||<^i(^)||- This essentially means that when the sectional 
curvature has a positive value then the radial geodesies tend to strongly converge/ focus, 
a fact which is expected/natural if one thinks in terms of General Relativity or any other 
generally covariant (diffeomorphism invariant) theory. Following the Rauch comparison 
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theorem it is easy to argue that a manifold with more negative sectional curvature will 
have geodesies diverging faster from each other than one with less negative, hence larger, 
sectional curvature. This is a different way of thinking about the geometric meaning of 
sectional curvature as compared to (3). The more difficult point is to prove that the vari- 
ations in curvature varying continuously from point to point and in different directions 
leave the rates of expansion and contraction of the appropriate stable and unstable dis- 
tributions essentially unchanged. To this end, one can use the rigidity of the hyperbolic 
space (manifold of constant negative curvature K < 0) and the properties of its ideal 
boundary [20], [29]. 

Now we can see why the largest Lyapunov exponent of a system described by the Tsallis 
entropy is zero. As was pointed out in previous work [30], [31] an effective framework in 
which the Tsallis entropy can be understood, and where nice algebraic properties with 
respect to the Tsallis entropy composition arise, is provided by constructing a deformation 
of the reals by the entropic parameter q which results in the set Rg. The field 
isomorphism : R — >■ is given by [30], [31] 

and we provided some preliminary comments on the induced metric properties of such Tg 
for q e (0, 1). We saw in [30] (see also [31]) that the importance of such a map may 
not so much its exact form, as much as that that it is of exponential form. What Tg 
does is that it maps distances on IR" to exponentials of such distances in M^. Prom the 
viewpoint of the Cartesian product all distances of R" appear to have shrunk, as 
they appear to be measured in a logarithmic scale. Equivalently, what Tg does is that 
it uses exponentially large scales to measure distances on any Euclidean space R", and 
by extension on any Riemannian manifold. Such an exponential scale which is used to 
measure distances is natural, from the viewpoint of the composition properties of theTsal- 
lis entropy. Then any exponentially increasing quantity, such as the separation of nearby 
geodesies on manifolds with negative sectional curvature, will appear to have a polynomial 
(power-law) asymptotic behavior. As a result, from the viewpoint of the metric struc- 
ture induced by the Tsallis entropy, all geodesies on Riemannian manifolds of negative 
sectional curvature separate with distances which increase at polynomial rates. Hence 
the resulting largest Lyapunov exponent of the geodesic flow is zero, a behavior coined 
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"weak chaos" [2] (and references therein). This is the kind of dynamical behavior that 
the Tsallis entropy is conjectured to encode. For preliminary indications that this may 
indeed be so, see [2], [11] and references therein. To our knowledge, no exceptions are 
currently known to the fact that the Tsallis entropy describes systems with zero largest 
Lyapunov exponent. Following the picture presented above, if someone wished to have 
positive largest Lyapunov exponent with respect to the metric induced by Tg from the 
Tsallis entropy, then nearby geodesies would have to separate at rates which are increasing 
at least as fast as e^^'°^*\ This cannot occur for manifolds of constant negative sectional 
curvature, as is immediately obvious from (8). 



4. Generalized geometric structures and discussion 

The preceding section justifies why the largest Lyapunov exponent is zero, for systems 
described by the Tsallis entropy. Hence if one persists in quantifying the asymptotic 
behavior of geodesies, then a modification of the definition of Lyapunov exponents is 
required. Since the exponential growth of nearby perturbations is a very strong condition 
as seen above, a polynomial (power-law) rate of increment may be more appropriate. To 
capture this rate of growth, modified Lyapunov exponents can be defined by 

logA,(y)= limhog\\dJtY\\ (12) 

Such modified exponents encode perturbations growing polynomially as t'^''^^\ This is 
in agreement with the conjectured behavior [2], [7]-[10] 

| = A,e (13) 

for systems described by the Tsallis entropy. Naturally, these modified Lyapunov expo- 
nents will explicitly depend on the value of the non-entropic parameter q. Since, as is 
conjectured, the systems described by the Tsallis entropy may possess more than just one 
value of the non-entropic parameter q [2], [11], depending on which property of the system 
is described, such modified Lyapunov exponents would explicitly depend on the value of 
q characterising the sensitivity of the system to infinitesimal perturbations ^sen [2] ■ 

For this kind of definition to be of physical interest though, there must exist physical 
systems exhibiting such a behavior. To be more precise, we would like to know whether 
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there are geometric structures which are the configuration or phase spaces of systems 
the stabihty of the geodesic fiow of which is characterized by the modified Lyapunov 
exponents. Such structures cannot be Riemannian manifolds, as was established in the 
previous section, since what we are seeking are systems whose perturbations grow accord- 
ing to a polynomial (power) law. Additionally, we would like to understand how do such 
structures arise in cases of physical interest. 

First, a generalization of the concept of minimal geodesic: let (M, d) be a metric 
space and let x, y E M. Then a geodesic between x and y is an isometric embedding 
c of the interval [0,d{x,y)] into M such that c(0) = x and c{d{x,y)) = y. To 
proceed, one makes simple but fundamental observation: when Riemannian manifolds 
have positive sectional curvature then triangles whose sides are locally minimal geodesies 
are "fatter" than corresponding triangles in Euclidean space. The triangles in both spaces 
are assumed to have corresponding sides of the same lengths. When such manifolds have 
negative sectional curvature then the corresponding triangles are "thinner" than in the 
Euclidean case. The comparison can be equally well be made with a space form of sec- 
tional curvature k instead with the Euclidean space [17] - [20]. When, for instance, 
a manifold has sectional curvature smaller than A; < (which is the case of interest in 
what follows), then its geodesic triangles will be "thinner" than those of the space form 
of sectional curvature k, when both triangles have the same corresponding side-lengths. 
A.D. Alexandrov [32], [33] reverses this argument by postulating the "fatter" or "thinner" 
triangles condition in order to define curvature in a comparison sense. What one gains by 
adopting such a definition is that it applies without a pre-requisite of smoothness of the 
underlying structure, nor of even a local Euclidean behavior such that of a topological 
manifold. A metric tree, namely a graph without loops, where each edge connecting two 
vertices is assigned a unit length is a standard example of such a structure, which does 
not look like a manifold at all, but for which curvature in the comparison sense can still 
be defined [34], [20]. What we lose by using such a definition is that we cannot speak 
about a space whose sectional curvature is of a particular value but only about a space 
whose sectional curvature is smaller (or larger) than a particular, negative (for the case 
of our interest), value. 

By "spaces" we will mean metric spaces with intrinsic length structures, also called 
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geodesic spaces, although more general constructions are possible [20]. The geodesic 
space requirement is that for any two points of such a space there is at least one mini- 
mal geodesic joining them. For the case of Riemannian manifolds of negative sectional 
curvature such a geodesic is actually unique [17]-[20], [24], [29]. Spaces with an upper 
bound G M on their curvature, which is defined in a comparison sense, are called 
CAT{k) spaces [34], [20]. In the case of metric trees, mentioned above, one can easily 
check that they are indeed CAT{0) spaces [34]. Other examples of CAT{k) spaces, 
are of course, Riemannian manifolds of sectional curvature bounded above by k. There 
is still some "rigidity" in CAT{k) spaces, else not many general statements could be 
made about them. Consider, for instance, a linear space endowed with a norm, which 
is also CAT{k) for A; G M. Then this turns out to be an inner product space. On 
the other hand, some properties of CAT{k) spaces may be unexpected. The local to 
global principle does not apply for such spaces, for instance: an infinitesimal curvature 
bound does not necessarily imply a global curvature bound. For a CAT{k), k < space, 
using the triangle inequality, one can prove that the asymptotic separation of geodesies 
is either bounded or exponential [33], [34]. Nothing new or unexpected occurs here when 
compared to the Riemannian case. But, in stark contrast to the Riemannian case, it is 
possible to have a quadratic deviation of geodesies in a CAT{0) space. The exact state- 
ments and the associated construction can be found in [351*1 In such cases the modified 
Lyapunov exponents defined in the previous section may be of interest in capturing the 
asymptotic behavior of the geodesic flow. In [36] we presented in more detail some of 
the above constructions and argued in greater detail, why CAT{k) spaces are indeed 
useful for describing interacting systems having different values of the entropic parameter. 

Why would such generalized metric structures be of interest to Statistical Mechanics? 
Shouldn't Riemannian manifolds and their asymptotic geodesic behavior be sufficient? 
The answer is two- fold: First, it is indeed true that configuration and phase spaces are 
very well modelled by Riemannian manifolds and their (iterated) tangent and cotangent 
spaces as well as their tensor products. A chance for something more general to be needed 
comes from at least two sources: one is when the symmetries of the system impose con- 
straints which introduce singularities to these Riemannian manifolds. A typical case is 
the non-free action of a group on the configuration/phase space. This is not an unusual 

*We would like to thank Professor Panos Papasoglou for bringing this reference to our attention 
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occurrence in multi-particle systems and is even more pervasive in the statistical descrip- 
tion of classical and quantum fields. Another case is the case of lattice models such as 
the Ising, Potts etc. There, the configuration space is discrete so it requires an approach 
somewhat different from the one applied for continuous systems. Another potential rea- 
son for use of genelarized geometric structures is the thermodynamic limit. The exact 
way of even appropriately defining the thermodynamic limit tends to be quite non-trivial, 
especially for systems with long-range interactions [2], [27]. It is not obvious that the 
thermodynamic limit of such a system will still be a Riemannian space. We can illustrate 
this with a simple, pictorial, example of a limiting process. Consider a two dimensional 
cone whose apex has been rounded. This is a Riemannian manifold. Consider now a se- 
quence in which a neighborhood of the apex gets more and more sharp. The limit of such 
a sequence will be the cone, which however is not a Riemannian manifold at its apex. 
Notice that we have avoided providing any details regarding the nature and definition 
of such a limit, by just alluding to a picture. Such an argument can be made though 
for a variety of limiting procedures, although the one that we have in mind in this case 
is the pointed Gromov-Hausdorff convergence [20], [37]. By analogy, we may be forced 
to consider generalized metric structures as the (thermodynamic) limits of Riemannian 
configuration/phase spaces. The second source of non- Riemannian geometric structures 
is an obvious one: the majority of dynamical systems for which a statistical description 
may be of interest are not Hamiltonian or even dissipative modifications of Hamiltonian 
systems. For these, the Tsallis entropy seems to be quite useful in encoding the asymp- 
totic (long-time) features of their orbit structure [2] even though such systems may not be 
described by smooth structures on Riemannian manifolds. Someone can readily mention 
the fractal-fike behavior in the phase space of systems fike these, which also motivated 
the introduction of the TsaUis entropy in the first place [1],[2]. After all, many (if not 
most) configuration spaces of systems of physical interest are actually spaces of mappings 
where it is the highly irregular mappings that provide the emerging behavior. An example 
of this case are the normal (or generalized) random walks, and by direct extension, the 
field theoretical treatment of quantization/ thermalization of systems via the functional 
integral formalism. 

A question directly related to the present work, is why the Tsallis entropy induces 
such an exponential increase in distances through which, in turn, results in the largest 
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Lyapunov exponent of the systems described by it to be zero. What is the root cause 
of such a behavior? Where exactly is this type of hyperboHc behavior, where distances 
are measured in an exponential scale, is encoded in the Tsallis entropy? Is this due to 
the particular metric properties of the field isomorphism that was constructed in 
[30], [31] or is this a generic behavior of the Tsallis entropy which just becomes manifest 
by such a construction of r^? The answer is the latter: the exponential expansion of 
the scale in which the distances are measured is already built in the Tsallis entropy. The 
construction of Tg [30], [31] is just one particular way to exhibit this and its associated 
properties. That this kind of hyperbolic behavior can be expected, is already encoded in 
the composition/generalized additivity property of the Tsallis entropy itself. How exactly 
this happens, a fact which is also loosely related to aspects of the uniqueness of the Tsallis 
entropy, will be discussed extensively in a near-future work. 

There is also a more general issue that may be of concern for some systems, such as 
systems with long-range interactions. In many such cases a cross-over behavior is observed 
in simulations of dynamical systems with few degrees of freedom [38] - [40] : such systems 
appear to be have power-law sensitivity to the initial conditions for a relatively long time, 
a situation dubbed as "quasi-stationary" . Such a behaviour can be described by the Tsal- 
lis entropy as ia well-known. After a long time passes though, such systems exhibit an 
exponential sensitivity the initial conditions, a behaviour that can be described by the 
BGS entropy. In essence one observes that such systems transit from weakly chaotic to 
strongly chaotic regions, depending on the elapsed time. Such cross-over in the sensitivity 
behaviour of the initial conditions begs for an explanation. The Lyapunov exponents are 
ill-adapted to explaining such a complex behaviour because they are asymptotic invari- 
ants. Hence their calculation should ignore any transient behaviour that may be present 
during the evolution of the dynamical system, no matter how long this finite transient 
behaviour may last. To deal with such a complex behaviour one may have to accept that 
the nonextesive sensitivity g^gn may be a function of the evolution parameter ( "time" ) . 
In the simplest case, and in order to retain some predictive power of this formalism, g^en 
should be chosen to be piecewise constant in each of the different temporal regimes of 
interest. It then becomes part of the theory to not only predict such values for qgen but 
also to determine the corresponding time intervals for which such values are effectively 
constant. It is unclear to us, at the present time, how much one would be able to proceed 
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along these lines and to what extent this approach would give results consistent with the 
numerical data, so it may be a direction of investigation worth pursuing. 

Acknowledgement: The author is grateful to Professor C. Tsallis for emphasizing 
that it is the largest Lyapunov exponent that vanishes in systems described by the Tsallis 
entropy, but not all of them, and for bringing to his attention reference [31] of which the 
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